Oscillating phantom in F(R) gravity 
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Abstract 

We investigate the oscillating effective equation of state (EoS) of the universe around the phantom 
divide in the framework of F(R) gravity. We illustrate the behavior of F(R) with realizing multiple 
crossings of the phantom divide. 
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Recent observations have supported that the current expansion of the universe is accel- 
erating 

HQ. 

There are two broad categories to account for this scenario 
One is the introduction of "dark energy" in the framework of general relativity. The other 
is the study of a modified gravitational theory, e.g., F(R) gravity, in which the action is 
represented by an arbitrary function F(R) of the scalar curvature R (for reviews, see j^, 7|). 
It is also known that, according to the observational data 8J, the ratio of the effective 
pressure to the effective energy density of the universe, i.e., the effective equation of state 
(EoS) w e ff = Pes/Pes, may evolve from larger than —1 (non-phantom phase) to less than —1 
(phantom phase 9J). This means that it crosses —1 (the phantom divide) at the present 
time or in near future. There are many models to realize the crossing of the pha ntom divide 



id, 



ii 



121 ]. an explicit 



(for a detailed review, see 4j). In the framework of F(R) gravity 6 
model with realizing a crossing of the phantom divide has been constructed in Ref. [13 1 

n 

and its thermodynamics has been examined [14|]. In this model, only one crossing of the 



phantom divide can occur because there appears the Big Rip singularity [15( at the end of 
the phantom phase. However, in the framework of general relativity multiple crossings of 
the phantom divide can be realized, e.g., in an oscillating quintom model [16| or a quintom 



with two scalar fields 



In this paper, we study the oscillating effective EoS of the universe around the phantom 
divide in the framework of F(R) gravity. We investigate the behavior of F(R) with realizing 
multiple crossings of the phantom divide. We use units of k-g = c = h = 1 and denote 
the gravitational constant 8irG by k 2 = 8n/Mp 2 with the Planck mass of Mpi = G^ 1 ! 2 = 
1.2 x 10 19 GeV. 

We start the reconstruction method of F(R) gravity proposed in Ref. 
F(R) gravity with matter is as follows: 

'F(R) 



. The action of 
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matter 
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where g is the determinant of the metric tensor g^ u and /^matter is the matter Lagrangian. 
By using proper functions P(4>) and Q(<p) of a scalar field 0, the action in Eq. (TjQ) can be 
rewritten to 



S= J d 4 Xv^|^2 [^(0)^ + Q(0)]+ ^matter I ■ (2) 

The scalar field may be regarded as an auxiliary scalar field because it has no kinetic term. 



From Eq. ([T]), the equation of motion of is given by 

Substituting = (f)(R) into the action in Eq. (j2J) yields the expression of -F(-R) as 

F(R) = P(<f>(R))R + Q(<t>(R)). (4) 

From Eq. (j2J), the field equation of modified gravity is derived as 

l -g, v [P((f>)R + Q(0)] - R, U P(<P) - ^DP(0) + V„V„P(0) + k 2 ^) = , (5) 

where V M is the covariant derivative operator associated with g^ iu , □ = g^V M V '„ is the 
covariant d'Alembertian for a scalar field, and j^™ atter ) j s the contribution to the energy- 
momentum tensor form matter. 

We assume the flat Friedmann- Robert son- Walker (FRW) space-time with the metric, 

ds 2 = -dt 2 + a 2 {t)dx 2 , (6) 

where a(t) is the scale factor. In this background, the components of (/i, v) = (0, 0) and 
(ji,u) = = 1, • • • ,3) in Eq. P read 

-6P 2 P(0(t)) - Q(cj)(t)) - QH dP ^ f)) + 2k 2 p = , (7) 

2« + *H d -^l + (AH + 6H 2 ) P(0(t)) + Q( m + 2, 2 P = , (8) 

where H = a /a is the Hubble parameter with ' = d/dt and p and p are the sum of the 
energy density and pressure of matters with a constant EoS parameter it>j, respectively, with 
% being some component of matters. After eliminating Q(4>) from Eqs. (I7j) and (jHJ), we obtain 

«^« + 2 ™ )) + k2( , + p) . . (9) 

The scalar field may be taken as = £ if it is redefined properly. By representing a(t) as 

a(t) = aexp (g(t)) (10) 

in terms of a constant of a and a proper function of g(t) and using if = dg(<f))/ (d<fr), we 
rewrite Eq. ([9]) to be 

^ 2 pw _ dm dm 

d0 2 c/0 C*0 c/0 2 W 

+ k 2 (! + ">0 Pia^ 3(1+ " l) exp [-3(1 + «;<) ^(0)] = , (11) 



where pi is a constant. Moreover, from Eq. (J7J), we get 



-6 



dg{<j>) 



P{</>) - 6 



dg(cf>) dP(cf>) 



Ik 2 P,a- 3(1+UH) exp [-3 (1 + w t ) g(<f>)} . 



(12) 



We note that if we redefine the auxiliary scalar field <fi by <fi = with a proper function 
$ and define P(ip) = P($(^)) and <3(y?) = Q($(^)), the new action 



,5' 



F(R) 



2k 1 



matter 



F(R)=P(<p)R + Q(<p) 



(13) 
(14) 



is equivalent to the action in Eq. (j2J) because F(R) = F(R). Here, p> is the inverse function 
of $ and we can solve <p with respect to R as <p — <fi(R) = $ _1 (0(-R)) by using (f> = 4>{R). As 
a consequence, we have the choices in </> like a gauge symmetry and thus we can identify 
with time t, i.e., 4> = t, which can be interpreted as a gauge condition corresponding to the 
reparameterization of = 4>((p) 13) . Hence, if we have the relation t = t(R), in principle we 
can obtain the form of F(R) by solving Eq. (fTTj) with Eqs. (j3J) and (fT2"l) . 

We also mention that the phantom crossing cannot be described by a naive model of 
F(R) gravity. To demonstrate the crossing, F(R) needs to be a double-valued function, 
where the cut could correspond to U7 c fr = ~1- However, the crossing can be performed by 
the extension of F(R) gravity, whose action is given by P(<f))R + Q{4>)- 

We consider a cosmology at late times. To illustrate the behavior of F(R) with realizing 
multiple crossings of the phantom divide, we investigate the case in which the Hubble rate 
H(t) is expressed by an oscillating function, given by [jij 



H — Ho + Hi sin vt 



v 



— H n , 

2 p 



(15) 
(16) 



where H p = 2.13/i x 10" 42 GeV 



20( with h = 0.70 2l| is the present Hubble parameter. Here, 



we take H = aH\ and H\ to be constants with a > 1 so that H can always be positive and 
the universe can expand. By using Eqs. ([IB]) and ( fTBT) and t p ~ 1/Pp, the present time t p 
is expressed by t p = (37r/2) v~ x and i?i = H v j (a — 1). We note that the behavior of H in 
Eq. f[T5"j) is only acceptable at late times. 
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FIG. 1: Time evolution of u> e fr for a = 10 with t = vt. 



In the FRW background, the effective energy density and pressure of the universe are 
given by p c fr = 3H 2 /k 2 and p c n = ~ (%H + 3H 2 J /k 2 , respectively. The effective EoS 
w^cfr = Pcff/Pefr is defined as jfj] 

2if 



W e ff = -1 



3# 2 ' 



(17) 



which implies that a crossing of the phantom divide occurs when the sign of H changes. 
In case of Eq. ()15|) . multiple crossings of the phantom divide can be realized. This means 
that the phantom phase is transient like the model in Ref. 11] and hence there is no Big 
Rip singularity. Hereafter, we take <fi = t. We show the time evolution of w e ff in Fig. 1 
with t = vt. In all figures, we take a = 10. From Fig. 1, we see that at the present time 
t p = (37r/2) v" 1 , the universe enters the phantom phase w e ^ < — 1 from the non-phantom 
phase w c fi > — 1. Thus, a crossing of the phantom divide occurs. 

From Eqs. ([TO]) and (USD with H = aH u H = dg(t)/(dt) and R = 6 (h + 2# 2 ), we 
obtain 



m 

a(t) 



H\ ( at cos vt 



exp 



a 



« rr I 1 

h Hi [at cos vt 

- 1 \ v 



(18) 
(19) 
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R = 6H 2 



v 

~H X 



cos vt + 2 (a + sin vt)' 



(20) 



where we have taken a = exp [—a (a — 1)] so that the present value of the scale factor should 
be unity. 

We define X = cosvt and solve Eq. (I2"0"j) with respect to X. If a is much larger than 
unity, we can neglect the term proportional to sin 2 vt in Eq. (120]) and therefore obtain the 
approximate solutions 



X(R) 







16a 2 2 



-2a 2 + 



R 



R 



(21) 



where 
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Hi 
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(22) 



In what follows, we use the upper sign in Eq. (1211) . 

For simplicity, we consider the case in which there exists a matter with a constant EoS 
parameter w = p/p. In this case, Eqs. (II ip and (Tl2|) are rewritten to 



(Aa 2 f3 2 + i)X±a/3 {4f3X + 3) Vl - X 2 - p (l - X 2 ) 



dP{R) 



+ 



2Fii/ 



(1 + w)pexp < — 3 (1 + w) 



a 



a 



dR 

j3 (a arccos X — X) 



+ XP{R) 

= (23) 



and 
Q(R) 



-6p 2 (a ± Vl-X 2 Y P(R) - 36/3 2 (a ± Vl - X 2 ) (4a/3X ^ Vl - X 2 



2k 2 
-pexp 



-3(l + w) 



a 



a 



+ (a arccos X — X) 



dP(R) 
dR 

(24) 



respectively, where the upper signs correspond to smut = +\A — X 2 , while the lower ones 
smut = — Vl — X 2 . We will concentrate on the latter case, while the former one will be 
remarked at the end. Here, p corresponds to the present energy density of the matter. 
In particular, we use the present value of the cold dark matter with w = for p, i.e., 
p = 0.233p c [J, where p c = 3H 2 / (8ttG) = 3.97 x 10~ 47 GeV 4 is the critical energy density. 
From Eq. (j4j) and the first relation in Eq. (1221) . we have 



2k 2 



2k 



(25) 
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FIG. 2: P(R) and Q(R)/v 2 as functions of R for a = 10 and p = 0.233p c 



To examine F(R), we numerically solve Eqs. fl23|) - fl25|) . 

In Fig. 2, we depict P(R) and Q(R)jv 2 as functions of R. The range of i? is given 
by 0.22 < R < 0.75 by solving Eq. (123]) numerically, corresponding to —1 < X < 1 in 
Eq. (j2~TT) . Here, we have taken the initial conditions as P(R = 0.22) = 1.0 and dP(R = 

0. 22) /(dR) = so that at a smaller curvature, F(R)/ (2k 2 ) could contain the term Rj (2k 2 ), 

1. e., the ordinary Einstein-Hilbert action. By using Eq. (|2"5]) . we show the behavior of 
F(R)/ (2k 2 ) in Fig. 3. We note that the qualitative behavior of F(R)/ (2k 2 ) in Fig. 3 does 
not depend on the initial conditions and that the quantitative values of F(R)/ (2k 2 ) do 
not also depend on the initial conditions strongly. Furthermore, we illustrate w e ff(R) = 

2 



-1 - [2/ (3(3)} X{R)/ ( a - y 1 - X 2 (R)j in Fig. 4. The time evolution of R is given in 
Fig. 5. From Figs. 4 and 5, we see that multiple crossings of the phantom divide can be 
realized. We note that the results in the all figures are shown by dimensionless quantities. 
We also note that the lower sign in Eq. fl2T|) leads to 0.59 < R < 1.11 with the results being 
qualitatively the same as those obtained by the upper sign above. 

From Fig. 3, we see that F(R) increases in terms of R around the present curvature 
R(t = = 12 [2/ (37r)] 2 = 0.54. This behavior is reasonable because in the Hu-Sawicki 
model [22j of F(R) gravity, which passes the solar system tests, F(R) increases around the 
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present curvature. We mention that such a behavior is typical for a general class of viable 
F(R) gravities introduced in Ref. 23| to which the Hu-Sawicki model belongs. This class 
of modified gravities can satisfy the solar system tests and unify inflation with the l ate- 
time cosmic acceleration. As viable models of F(R) gravity, e.g., the models in Refs. 



24 



25 



26 



27 



28 



291 ] are also known (see also 30j). We remark that it is impossible to 



state anything on whether the reconstructed F(R) gravity can pass the solar system tests 
and cosmological constraints at early times because the reconstructed F(R) gravity cannot 
describe the behavior of F(R) gravity for large values of R. 

We note the stability for the obtained solutions of the phantom crossing under a quantum 
correction coming from the conformal anomaly. In Ref. it has been shown that the 
quantum correction of massless conformally-invariant fields could be small when the phantom 
crossing occurs and therefore the solutions of the phantom crossing could be stable under 
the quantum correction, although the quantum correction becomes important near the Big 
Rip singularity. In the present model of F(R) gravity with realizing multiple crossings of 
the phantom divide, the obtained solution can be stable because the phantom phases are 
transient and there is no Big Rip singularity. 



Finally, we would like to remark that for the other possible solution of sin vt = +yl — X 2 , 
i.e., the case of the upper signs in Eqs. (1231) and (124]) . F(R) approaches minus infinity as the 
curvature becomes larger than the present one although it increases in terms of R at smaller 
curvature than the present one. Such a behavior is incompatible with the Hu-Sawicki model, 
in which F(R) ~ R + const, at a much larger curvature than the present one. 

In summary, we have studied the oscillating effective EoS of the universe around the 
phantom divide in F(R) gravity. In particular, we have analyzed the behavior of F(R) 
gravity with realizing multiple crossings of the phantom divide. Our result can be interpreted 
as an explicit example to illustrate that multiple phantom crossings can occur in F(R) 
gravity as the scalar field theories such as an oscillating quintom model in the framework 
of general relativity. 
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FIG. 3: Behavior of F(R)/ (2k 2 ) as a function of R. Legend is the same as Fig. 2. 
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FIG. 4: Behavior of w e s (R). Legend is the same as Fig. 2. 
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FIG. 5: Time evolution of R. Legend is the same as Fig. 1. 
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